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ABSTRACT 

We present a Large Binocular Telescope (LBT) imaging study that characterizes the star cluster 
co mponent of nearby galaxy outer disks (beyond the optical radius i?2s)' Expanding on the pilot project 
of Herbert-Fort et al. (2009), we present deep (~ 27.5 mag U-band point-source limiting magnitude) U- 
and V-band imaging ofsixgalaxies: IC 4182, NGC 3351, NGC 4736, NGC 4826, NGC 5474, and NGC 
6503. We find that the outer disk of each galaxy is populated with marginally-resolved star clusters with 
masses ~ 10 3 M Q and ages up to ~ 1 Gyr (masses and ages are limited by the depth of our imaging 
and uncertainties are large given how photometry can be strongly affected by the presence or absence 
of a few stars in such low mass systems), and that they are typically found out to at least 2i? 2 5 but 
sometimes as far as 3 to 4 R25 - even beyond the apparent H I disk. The mean rate of cluster formation 
for li?25 < R < l-5i?25 is at least one every ~ 2.5 Myr and the clusters are spatially correlated with the 
H I, most strongly with higher density gas near the periphery of the optical disk and with lower density 
neutral gas at the H I disk periphery. We hypothesize that the clusters near the edge of the optical disk 
are formed in the extension of spiral structure from the inner disk and are a fairly consistent phenomenon 
and that the clusters formed at the periphery of the H I disk are the result of accretion episodes. 

Subject headings: galaxies: star clusters - galaxies: structure - methods: statistical - facility:LBT 



1. INTRODUCTION 

In the current paradigm of galaxy evolution, material 
continually accretes onto galaxies. As such, the outer 
extremities are particularly interesting environments, but 
they are notoriously difficult to study. The surface bright- 
ness of the diffuse stel lar component drop s well below the 
background sky l evel ( Pohlen et al.||2002 1 , the neutral gas 
becomes ionized (Maloney 1993), and observing the dif- 
fuse ionized gas requires long exposures on the largest tele- 
scopes while still b eing limited to within the ra dial extent 
of the neutral gas ( Christlein fc Zaritsky||2008 1 . 

The GALEX mission ( |Martin et al.||2005p highlighted 
an alternative approach to the study of outer disks by 
clearly identifying large populations of apparent young 
stell ar clusters in the outskirts of some nearby disk galax- 



ies dThilker et al.||2005| |Gil de Paz et~alp005] [Zaritsky fc 
Christlein|20'071 rCALEX UV imaging proved particularly 
useful in this regard because it provides the stark color 
contrast required to easily differentiate young, blue clus- 
ters from the sea of redder background objects. Of course, 
any method that distinguishes between young clusters and 
background objects can be used to the same effect, and so, 
narrowband Ha imag ing both predates the GALEX work 



(Fer guson et al.|1998[ ) and still provides new results (|Werk 
et al. 2010[ ). While having certain advantages, the draw- 
back of the Ha work, and to a lesser extent the GALEX 
studies, is that they focus on the youngest clusters, which 
limits the number of such objects available for study and 
also our ability to measure the long-term history of this 
galactic component. Nevertheless, the value of these stud- 
ies is evidenced by the increase in corresponding theoret- 



ical studies (e.g. Bush et al. 2008 Roskar et al. 2008a|b 



A 



2009) 



Kazantzidis et al 

As useful as GALEX has been for detecting large pop- 
ulations of outer disk clusters, this paper focuses instead 
on broadband optical imaging of outer disks for several 
reasons. First, the spatial resolution of GALEX imaging 
(~ 5" FWHM) is roughly six times poorer than what we 
normally achieve using ground-based optical telescopes. 
Therefore, the candidate clusters identified by GALEX 



are typic ally blends of multiple clusters (Gil de Paz et 



al. 2005), as can also be seen in the comparison we pro- 
vide in Figure [T] Until we establish what the typical mass 
of our detected objects is, we will refer to these objects 
as stellar knots, which is a purely observational definition 
of an individual object identified in the available images, 
to avoid prejudging their true physical nature. Of course, 
even after we have an estimate of the typical mass one 
must remain aware that any particular detection may be 
quite different from the mean. Second, our knot mass 
limit is an order of magnitude lower than that of the 
GALEX samples, leading to a corresponding increase in 
the number of identifications. Third, optical colors can, 
in principle, differentiate stellar populations that are sev- 
eral times older than those differentiated with the UV col- 
ors, thereby extending the baseline over which the phe- 
nomenon can be studied, and again increasing the sample 
size. If we can overcome the difficulty in distinguishing be- 
tween stellar knots and background sources in optical im- 
ages, high-quality ground-based data will provide a much 
greater number of knots, and potentially a larger sample 
of galaxies with such data, for any statistical analysis of 
outer disk properties. 



Fig. 1. — Comparison of GALEX (left) and LBT imaging (right) of a 3 arcmin-wide (4.1 kpc) region in the outskirts of NGC 4736. 



The key to success lies in identifying some property, in 
addition to color, that can be used to disentangle the stel- 
l ar clu sters from background galaxies. |Herbert-Fort et &L\ 
( 2009 ) used the self-clustering of the outer disk knots to 
that effect and presented results for NGC 3184. Here we 
apply both the "classical" color approach and our self- 
clustering method to characterize these outer disk popu- 
lations. With results drawn from these two approaches, 
we address questions regarding the ubiquity of outer disk 
knots, measure the knot formation rate, probe the radial 
extent of the knot distribution, and compare the distribu- 
tion of knots to the H I distribution. 

We present the results of a statistical study of nearby 
(< 15 Mpc) outer disks, using the 2 x 8.4m Large Binocular 
Telescope (LBT; Hill et al. 2006) and wide-field, prime- 
focus Large Binocular Cam eras (LBC; Ragazzoni et al. 



2006| |Giallongo et al.||2008[ ). We describe our data re 
duction, develop analysis tools, and apply these to deep 
£/-band and V-band imaging data of six galaxies. We 
demonstrate how the two separate statistical methods en- 
able us to trace outer disk cluster-like objects to large 
radii. In §2 we describe our observations, data reductions 
and source detections. In §3 we present color-magnitude 
diagrams (CMDs) of candidate outer disk sources, quan- 
tify the range of properties consistent with the observed 
knots, and characterize the knot populations around galax- 
ies. In §4 we describe the application of our restricted 
three-point correlation analysis of LBT knots, a similar 
analysis of GALEX knots, and a cross-correlation analysis 
of LBT knots and the underlying H I distribution. First, 
we present the different analyses, and then discuss them 
jointly on a galaxy-by-galaxy basis in §5. We summarize 
our results and conclusions in §6. 

2. OBSERVATIONS TO FINAL SOURCE CATALOGS 

We present new observations of six galaxies (IC 4182, 
NGC 3351, NGC 4736, NGC 4826, NGC 5474, NGC 6503) 
with the LBC-Blue imaging camera on the LBT on the 
dates listed in Table[l] These systems comprise a represen- 
tative sample of nearby disk galaxies by spanning a range 
of physical disk sizes. All are nearby (D ~ 15 Mpc) and 
only two, NGC 4826 and NGC 6503, are seen at high incli- 
nation (> 60°). For the four low inclination systems, their 
small distances and low inclinations allow us to measure 
the clustering of sources in the disk as a function of galactic 



radius. We construct the final galaxy mosaics from many 
individual 164-second dithered exposures through the U 
and V Bessel filters (hereafter U and V), except for NGC 
5474, for which we combine 33-second dithered exposures 
in U to avoid saturating a bright star in the field. The 
typical seeing during our observations was 0'.'8 and never 
exceeded - 1'.'3. IC 4182, NGC 4736, and NGC 5474 were 
observed under photometric conditions. We obtained sin- 
gle photometric exposures for NGC 3351, NGC 4826, and 
NGC 6503 on 2/8/08 (U k V), 4/20/07 (U & V), and 
4/22/07 (U) & 4/7/08 (V), respectively, for photometric 
calibration. We observed six photometric Landolt stan- 
dard star fields ( |Landolt 1992 ) for flux calibration and all 
magnitudes are on the Vega system. 

We follow the im age processing steps described by |Herbefo - 
Fort et al. (2009). We use a set of Interactive Data Lan- 
guage (IDIJM) scripts we created to construct 'master bias' 
frames that correct for global changes and the two dimen- 
sional structure in the bias level across the CCDs. We cor- 
rect for sensitivity variations across the CCDs by combin- 
ing dithered, twilight-sky flat-field images that are normal- 
ized by the four-chip median value of each bias-corrected 
flat field exposure. We then median-combine the normal- 
ized fiats to create the master flat frames and, finally, di- 
vide the bias-corrected science frames by the master flats 
to complete the processing of the individual science expo- 
sures. 

The target galaxies are well centered on the mosaic, so 
we sample the background of each exposure in twenty re- 
gions around the edges of the detector array. From these 
twenty regions we are able to obtain sufficient samples 
from which to judge the uniformity of the background and 
obtain and estimate of the uncertainty in our background 
estimate. The background is estimated using the IDL rou- 
tine MMivj^] which uses an iterative process to determine 
the mode of the sky values after rejecting outlying pixels 
(caused by stars, saturation, cosmic rays, hot pixels and 
bad columns). This step is done for each four-chip ex- 
posure and the minimum estimated sky from the various 
regions is subtracted from each corresponding chip. We 
then create a weightmap for every four-chip exposure to 

1 Develope d by Rese arch Systems, Inc. and owned by ITT; 
|http:/ /www.ittvis.com/ProductServices/IDL.aspx| 

2 Part of the Goddard IDL library, maintained by W. Landsman; 
http://idlastro.gsfc.nasa.gov/ 
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mask bad columns (detected by eye in the master flats), 
hot pixels and cosmic rays (detected using the IDL routine 
reject_CR, which finds features sharper than the PSF) 
when creating the mosaic images. 

To combine the individual frames into final exposures 



we first use SExtractor (Bertin & Arnouts [1996 ) to cre- 
ate a source catalog for every exposure so that SCAMP^] 
( |Bertin| [2006) can correct for optical distortions in LBC- 
Blue by solving the astrometry of the dithered exposures 
and creating smooth distortion maps for use b y SWarrf] 



In cases where the SDSS fifth data release (DR5; Adelman- 
McCarthy et al.||2007[ ) does not cover the scienc e field". 



we use the USNO-B catalog ( |Monet et aL||2003| ) as the 
reference. We use SWarp to mean-combine the spatially- 
aligned frames after accounting for the bad pixel maps. 
Our final mosaic images are rebinned to a ~ 0'.'22 pixel -1 
spatial scale and have their photometry normalized to a 
fixed counts arcsec -2 second -1 standard. The mosaics are 
flat to 0.5 — 1% and have total integration times of ~ 25 
minutes in the deepest areas (see Table [TJ. 

We build aperture photometry catalogs using SExtrac- 
tor, with source detection and aperture placement based 
on the deeper U-band mosaic. We select sources by iden- 
tifying groups of five or more pixels each with flux > la- 
above the background (so detections are in total > 3<r). 
Because we detect the integrated light from members of 
stellar groups, the 'knots', any photometric algorithm that 
requires a uniform obje ct shape for extraction is not op- 
timal for this work. See Herbert-Fort et al. (2009) for all 
of the chosen SExtractor parameters. We detect nearly all 
visually discernable objects beyond the optical radius R25, 
though our catalog becomes noticeably incomplete inside 
~ 0.8i?25- Any algorithm will have difficulty detecting 
sources over the bright, extended emission of the inner 
disk. Issues regarding the completeness relative to our vi- 
sual detections are negligible because the faint objects that 
are of interest in such a discussion are later rejected on the 
grounds of the low precision of their color measurement. 

We use processed photometric exposures of Landolt stan- 
dard star fields, taken on the same night as our individ- 
ual photometric exposures of the galaxies, for calibration. 
We account for an airmass and a color term when flux- 
calibrating on the standard Vega system. Finally, we boot- 
strap the photometry of the photometric exposures to the 
deep mosaics using ~ 10 stars common to both. 

Colors and magnitudes are measured using circular aper- 
tures. Colors quoted throughout are from apertures with a 
diameter fixed to four pixels (0'.'9, just larger than the typ- 
ical 0'.'8 FWHM of detected sources) , while V magnitudes 
are apertures with a fixed 10 pixel diameter (2724) that 
are then aperture corrected using stellar curves of growth 
(see Table [l] for those corrections, V acorr ). Aperture cor- 
rections were calculated from ~ 8 isolated, unsaturated 
stars measured in 15 apertures spanning 2 — 50 pixels in 
diameter (or 0'.'4 — 11'! 2). We set SExtractor to mask and 
correct for contaminants. We check our photometry by 
comparing the U and V apparent magnitudes of ten well- 
isolated objects across the fields with those provided by 
SDSS-DR5 when available, otherwise USNO-B, converted 
from u, g, and r to either U or V using the transformations 

3 Version 1.4.0; [ http: //terapix. iap.fr/soft /scamp 



of Jester et al. (2005). Aside from the systematic offset in 
our V magnitudes from NGC 4826 (see below) , our results 
are consistent with the transformed SDSS photometry to 
within the transformed photometric errors. 

The night of April 24, 2007 had poorer seeing than usual 
(~ l'/3) during our V exposures of NGC 4826. As a result, 
the U — V colors of NGC 4826 sources, as estimated by eye 
from the CMDs presented in the next section, appear to 
be artificially blue by ~ 0.5 mag relative to the source dis- 
tributions of the other galaxies. This shift is as expected 
given the large aperture correction in V (Table fl]) and the 
fact that the colors are aperture-matched (and nence not 
aperture corrected). A similar color offset, though of lower 
magnitude, is also just noticeable in the color-magnitude 
diagrams of NGC 5474 and NGC 6503 and these galaxies 
have the next highest aperture corrections. The solution 
to this problem is to PSF-match all the frames, but that 
would result in a significant degradation of the data. We 
choose not to apply the PSF matching because these color 
offsets do not impact our results. None of the subsequent 
discussion is predicated on precisely measured colors and 
NGC 4826 (the galaxy showing the strongest effect), as 
well as NGC 6503, are eventually dropped from the anal- 
ysis because of their relatively high inclinations [i = 61° 
and 74°, respectively). 

To complete our final source catalogs, we reject sources 
whose internal U — V color error is > 0.5 mag (magni- 
tude errors are provided by SExtractor and propagated 
in the standard manner) and interactively mask regions 
around bright stars (those showing diffraction spikes and 
halos of scattered light) because the SExtractor catalogs 
have artificially low counts in those regions, which creates 
artificial structure in our correlation maps. In summary, 
our catalog is surface brightness limited due to the crite- 
rion on flux per pixel, magnitude limited due to the cri- 
terion on number of pixels above the particular surface 
brightness, and color limited due to the requirement of a 
modest uncertainty in the color measurement. The latter 
two in particular drive some of the structure seen in the 
color-magnitude diagrams discussed in §3. Our catalog is 
also subject to the confusion limit, although evidently at a 
finer resolution limit than GALEX (see Figure]!]). We pro- 
vide the number of sources in our final catalogs between 
1.0 - 1.5i? 25 , 1-5 - 2.0i? 25 , and 2.0 - 2.5i? 2 5 and the area 
in each annulus after masking (Table 2]). The chosen outer 
limit in this Table, 2.5i?25, is somewhat arbitrary, loosely 
based on previous o uter disk studies, for example Za ritsky| 
& Christleinl (120071. 



4 Version 2.17.1; |http://terapix.iap.fr/soft/swarp 



3. OUTER DISK CLUSTER POPULATIONS 

We will eventually appeal to our clustering analysis to 
tease out the most information regarding the distribution 
of the knots in the outskirts of these galaxies. However, it 
is worthwhile first examining color-magnitude diagrams to 
gain some intuition regarding the nature of the knots. The 
aim of this section is therefore not a detailed description 
of the knots, because we are dominated by uncertainties in 
accounting for the dominant background population, but 
rather a broad appraisal of the knots. Estimated masses 
and ages are highly uncertain. 

In Figure [2j we present the color-magnitude diagram, 
CMD, of sources between 1.0 — 1.5i?25 in the field of IC 



4182 and tracks from default Starburst99 (Lcithercr et al. 



4 



1999 |Vazquez fc Leitherer 2005 ) that show the locations of 
model clusters of tixed mass and solar metallicity for ages 
between 1 Myr to 3 Gyr. We have also examined models 
with l/5th and l/50th solar metalliciites. The l/5th solar 
model is indistinguishable at the level of precision relevant 
for the qualitative conclusions we reach (< 0.5 mag), and 
the 1/50 solar model differs significantly, > 1 mag, but 
only for ages < 10 Myr, at colors where we find few knots 
in most of our galaxi es. Given current observational con- 



strains (Werk et al. 2011), we do not expect outer-disk 
clusters to have extremely sub-solar metallicities. The up- 
per and lower tracks represent 10 4 M Q and 10 2 M Q clus- 
ters, respectively, scaled from a simulated 10 6 M Q cluster 
that adequately samples the upper mass range of a Kroupa 
IMF. Because we ignore t he stochastic sampling of the 
IMF at low cluste r masses ( Cervino fc Luridiana|2004 Fa- 
gioli ni et al.|2007| ) , a simple comparison between data and 
the scaled model tracks can lead to an underestimate of 
cluster mass and an overestimate of the age. Therefore, 
these tracks are only meant to provide a general impres- 
sion of the cluster masses and ages consistent with the 
range of sources in our sample. 

More instructive are the background-subtracted Hess di- 
agra ms (Figures 3 — 8) that we create in the same manner 
as in 



Herbert-Fort et al. 



(2009), to statistically constrain 
the color and magnitude range of sources most likely to 
be clusters associated with each galaxy. A Hess diagram 
is a plot of the number of sources within color-magnitude 
bins (we use square bins of 0.2 mag here). For these par- 
ticular diagrams, we create both a 'background' Hess di- 
agram from a region far outside the galaxy and a 'outer 
disk+background' Hess diagram from the region of inter- 
est (say 1.0— 1.5i?25), scale the counts in the 'background' 
Hess diagram by the relative areas of the two regions, and 
subtract it, bin-by-bin from the 'outer disk+background' 
Hess diagram. The three panels in each Figure show the 
residual source density between 1.0 — 1.5 R25 (left; dark 
regions are positive counts), 1.5 — 2.0 R25 (middle), and 




Fig. 2. — CMD of all sources surrounding IC 4182 (including back- 
ground sources) between 1.0 — 1.5i?25 (after masking areas around 
bright stars; see text), with median lcr errors as a function of ap- 
parent magnitude shown on the left. A reddening vector corre- 
sponding to 1 magnitude of extinct ion in V, calculated using re- 
sults from Riekc & Lebofsky (1985), is shown at lower right. The 
work oflAlberts et al.|12Ulip suggests low extinction in outer disks 
(E(B—V ) < U.3). The tracks represent solar- metallicity 10 3 (lower) 
and 1O 6 M0 (upper) star clusters as a function of age. These tracks 
do not account for stochastic sampling of the stellar mass function. 
See §3 for more details. 



2.0 — 2.5 i?25 (right). Solid and dotted black contours out- 
line signal above and below the background, respectively, 
at > 90% confidence level (CL), calculated using the low- 
count, Poisson single-sided upper and lower limits from 
Gehrels (1986). Because the solid contours surround pix- 
els whose individual value is above the background at the 
90% CL, groups of such pixels are detections at a much 
greater CL. Overplotted is a 10 3 M Q Starburst99 model 
cluster, scaled from a 10 6 M e cluster, as shown in Fig- 
ure [2] See Table [2] for source counts and other details of 
the Hess diagrams, including our estimates of the effective 
surface brightness of the outer disk cluster components. 

All six galaxies show excesses in their background sub- 
tracted Hess diagrams between 1.0 — 1.5 R25, roughly trac- 
ing the 1O 3 M cluster track. The tilted contours in the 
lower portion of all the panels reflects the color selec- 
tion (bluer knots are visible to fainter V magnitudes). 
All galaxies except NGC 5474 show suggestive excess be- 
tween 1.5 — 2.0 R25, though the noise is noticeably higher. 
We are skeptical, although not dismissive, of apparent ex- 
cess between 2.0 — 2.5 R25 because those diagrams are so 
strongly peppered with oversubtraction (nevertheless, the 
morphology of signal beyond R25 does follow that seen 
within i?25)- The principal source of uncertainty here are 
the variations in the background population. The effect of 
this can be seen in the appearance and disappearance of 
regions of oversubtraction, and in the change in location of 
such oversubtraction, within the various Hess diagrams. In 
contrast, the excesses seen in all six galaxies straddle the 
model track, lending credence to our association of such a 
signal with a physical, associated population of sources. 

In Figure [9] we present the average Hess diagram for our 
six galaxies, with populations combined at constant My, 
with and without NGC 4736. We single out NGC 4736 
because it is the nearest and largest galaxy (in angular ex- 
tent) in our sample, resulting in a disproportionate num- 
ber of detected sources and a less well-determined back- 
ground level. These average Hess diagrams, either with 
or without NGC 4736, show convincing excess out to at 
least 2.0 R25 in that there is more area enclosed within the 
solid contours than within the dashed contours (alterna- 
tively, based on Table [2] each galaxy has an average excess 
of over 70 knots between 1.5 and 2 R25, even if we exclude 
NGC 4736). We conclude that there is indeed a population 
of knots extending well beyond R25 in all of these galax- 
ies and that the knots are consistent with a population 
of stellar clusters. We will now refer to these objects as 
stellar clusters. The exact distribution of masses and ages 
is difficult to disentangle from these plots given the uncer- 
tainties in the luminosities and colors of low mass clusters 
arising from the stochastic sampling of the stellar mass 
function, the uncertainties in the modeling, and the large 
uncertainties arising from the background subtraction. 

3.1. An Estimate of the Cluster Formation Rate 

To gain some intuition on the implied cluster formation 
rate from these diagrams and to further test our assertion 
that these are clusters, we simulate a population of 1O 3 M0 
clusters forming at a specified constant rate over the pre- 
vious several Gyr and plot the resulting Hess diagrams, 
accounting for the photometric uncertainties as estimated 
from the data for NGC 4736 (but not the uncertainties in 
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Fig. 3. — Background-subtracted Hess diagrams for IC 4182, made from sources between 1.0 — 1.5 R25 (left), 1.5 — 2.0 R25 (middle), 
and 2.0 — 2.5 R25 (right). Dark regions represent positive counts. Solid and dotted black contours outline signal lying above and below 
the background at a > 90% confidence level (CL), respectively; the dotted contours show any oversubtraction. Ovcrplotted is a 
Starburst99 model cluster, scaled down from a 10 6 A/q cluster. Symbols represent the age of the cluster along the model sequence. 
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Fig. 4.— Same as Figure 3 but for NGC 3351. 
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Fig. 6.— Same as Figure 3 but for NGC 4826. 
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Fig. 7.— Same as Figure 3 but for NGC 5474. 
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Fig. 8.— Same as Figure 3 but for NGC 6503. 
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Fig. 9. — The average Hess diagram for all six galaxies, with sources combined as a function of My, including and excluding NGC 4736 
(top and bottom panels, respectively). The three radial ranges are the same as in the individual Hess diagrams (1.0 — 1.5 R25, 1-5 — 2.0 R25, 
and 2.0 — 2.5 i?2s). Dark regions are positive counts. Solid and dotted black contours outline signal lying above and below the background 
at a > 90% confidence level (CL), respectively; the dotted contours show any ovcrsubtraction. Overplotted is a 10 3 A4q Starburst99 model 
cluster, scaled down from a 10 6 Mq cluster. There is excess signal out to at least 2 /?25- 
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Fig. 10. — Simulated Hess diagrams made from a Starburst99 10 3 Mq cluster track, with CMD uncertainties corresponding to those found 
for NGC 4736, for different cluster formation rates (top left: one every 10 X 10 6 yr; top right: one every 10 6 yr; bottom left: one every 10 5 
yr; bottom right: one every 10 4 ). No background is considered when making the simulated diagrams. 
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the background subtraction nor the stochastic sampling of 
the stellar mass function, in a figure analogous to Figure [9] 
for four different cluster formation rates. We assume no 
cluster dissolution or disruption, so this comparison will 
provide a lower limit on the cluster formation rate. Given 
the uncertainties in applying these models to the data, 
arising from the causes outlined above, this calculation is 
intended only as a plausibility argument for relating the 
knots to stellar clusters. 

Comparing Figures [9] and 10, we conclude that the data 
exclude cluster formation rates significantly higher than 
one every ~ 10 6 years for clusters of mass > 10 3 M Q (if 
there is no significant cluster dissolution) because we do 
not detect a population of blue sources at My < —7 
mag. This conclusion is relatively insensitive to the is- 
sue of stochastic sampling because that becomes less of a 
factor for the more massive clusters. On the other hand, 
rates significantly lower than 10~ 6 yr _1 fail to produce a 
sufficiently significant population of sources. This range of 
rates is consistent with a calculation where we simply take 
the number of outer disk knots in the 1.0 — 1.5 R25 Hess 
diagrams (~ 400), the maximum age of a 1O 3 M cluster 
in our diagrams (~ 1 Gyr), and a uniform rate of forma- 
tion over that time, which results in a rate estimate of 
one cluster every ~ 2.5 Myr (4 x 10~ 7 yr _1 ). Converting 
this cluster formation rate to a stellar cluster mass rate 
implies that ~ O.OO4M pc~ 2 Gyr -1 is being tied up in 
stellar clusters at these radii (assuming i? 2 5 = 5 kpc). 

Alternative estimates of t he outer disk cluster for mation 
rate exist for comparison. |Ferguson et al.| ( |1998[ ), using 
deep Ha imaging, measured the outer disk star formation 
rate densities of NGC 628, NGC 1058 and NGC 6946 to 
be between ~ 0.01 - O.O5M pc" 2 Gyr" 1 . There are at 
least four potential explanations for our significantly (lOx ) 
lower formation rate, all of which probably contribute to 
the difference: 1) the cluster format ion rates are highly 
stochastic and Ferguson et al. (1998) happened to catch 
these galaxies in an elevated phase relative to ours, 2) clus- 
ters disassociate, and so many are missing in our sample 
which samples older clusters, and 3) the Ha technique, 
which is sensitive to low mass clusters, is measuring a dif- 
ferent mass component that happens to contain a larger 
fraction of the total mass, and 4) our estimates of the 
masses are significantly corrupted by our neglect of the 
stochastic effects and modeling uncertainties. 

Regarding the first possibility, we know from a com- 



parat ive study of GALEX knots (Zaritsky & Christlein 
2007) that only a fraction of galaxies (~ 25%) show sig- 
nificant overdensities of bright blue knots. If all galaxies 
have an outer disk population, as we seem to find on the 
basis of optical imaging pre sented here and spectroscopy 
( |Christlein fc Zaritsky|200"8 ), then the formation rate must 
be highly variable (with a duty cycle of about 25% for 
GALEX-detectable knots). The effects of this stochas- 
ticity should be even more dramatic in Ha. Regarding 
the second option, we know that in certain environments 
where we have clusters spanning a range of ages, and can 
therefore do the study, that only a small fraction of all 
star clusters survive and we believe we u nderstand the 
physical mechanism for this evolution (see Spitzer 1958 
Gieles et al. 2011 and many related studiesj! Depend 



tidal stresses) , the rate of cluster dissolution may be lower 
in the outer disks, but unlikely to be negligible. Davidgc 
et al. ( |2011 ) find that clusters in the outer disk of M33 
dissipate on a time scale of 100 Myr. The clusters we de- 
tect are ~100 Myr old and older, and so likely to be a 
remnant population. This possibility is g i ven fu rther sup 
port by the recent work of 
set of five galaxies in whic 



Alberts et al. (2011), who in a 
they are able to measure an 
age distribution of a set of massive (M ~ 10 5 — 1O 6 M0), 
outer disk clusters, find that the age distributions are all 
peaked toward early times (~ 100 Myr and often within 
their innermost age bin of 50 Myr) even though they sam- 
ple to ages of 1 Gyr . Finally, regarding the third option, 
Davidge et al. (2011 ) calculate that most of the outer disk 
clusters in M33 form at lower masses (50 — 250M Q ) and so 
it may be the case that many of the Ha-detected clusters 
are of similar masses and below our detection threshold. 
It may also be the case that many of our clusters are low 
mass and boosted by stochastic effects into detectability. 
As such, any quantitative determination of the mass func- 
tion of these clusters will require simulations that include 
detailed treatments of such effects, as well as of dynami- 
cal evolution and selection. Importantly, if we are missing 
90% to 95% of the outer-disk clus ter formation, as sug- 



mg on tEe driver for cluster dissolution (mass loss versus 



gested by the Ferguson et al. ( 1998 ) results, then we must 
bear in mind that all of our subsequent estimates of stellar 
mass densities at these radii need to be multiplied by a fac- 
tor of 10 to 20. For now, we ascribe the discrepancy to one 
of the effects described above rather than to a wholesale 
missing population from our catalog. 

If the outer disks in our sample account for a star forma- 
tion rates of 4 x 10~ 4 M Q yr _1 (one 1O 3 M0 cluster every 
2.5 Myr for 1.0R 25 < R < 1.5R 25 ), then - 4 x 1O 6 M of 
stars have formed in that annulus over the lifetime of the 
galaxy (taken as 10 10 yr). As we discussed above, this is 
a lower limit since it ignores cluster dissolution and po- 
tential selection effects. If we multiply this number by 20, 
then we conclude that this limited region of the outer disk 
could contain as much as 10 8 M© of stars, or about 1% of 
the stars in a typical large spiral. 

We conclude that adopting rough estimates of the typi- 
cal mass and age of these clusters results in an outer disk 
star formation rate that is consistent with other estimates 
in that it lies below those other estimates. As such, our 
claim that these are indeed outer disk star clusters does 
not conflict with those other observations. 

4. CLUSTERING IN THE OUTER DISKS 

The analysis of the color-magnitude diagrams are lim- 
ited by the immense contamination from background ob- 
jects. As we mentioned before, while at UV wavelengths 
that contamination is held in check, at optical wavelengths 
another method is needed to help differentiate outer disk 
clusters from background objects. To enhance the contrast 
between outer disk clusters and the background, we now 
utilize spatial correlations among the clusters and between 
the clusters and other detected, outer-disk components. 
First, we describe the self-correlation of LBT-knots. We 
will discuss the results for individual galaxies in detail, 
and the sample as a whole, in §5. Second, we present 
the application of the same technique to existing GALEX 
data for these same galaxies. In general, there are fewer 
GALEX knots per galaxy, so the statistical information is 
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poorer, but those data offer an independent check on our 
LBT results. Finally, we cross-correlate the position of the 
LBT-knots with the H I in an attempt to confirm tenuous 
evidence for very distant clusters through their association 
with neutral hydrogen at these large radii. 

4.1. Self-clustering of LBT Knots 



Following Herbert-Fort et al. ( 2009 1 , we present restricted 
three-point correlation maps to trace the self-clustering of 
knots in the outer disks (see Figure 11 for a description of 
the radii used). The self-clustering of outer disk clusters 
provides an enhanced contrast relative to the background 
(which has a different angular correlation function). In- 
stead of measuring a radial profile of detected sources and 
subtracting some average background level, which is sim- 
ply related to the a zimuthally-averaged two-po int corre- 
lation function (e.g. |Zaritsky fc Christlein]|2007 [), here we 
use the self-clustering of knots to highlight regions with 
clusters and remove signal from large-scale background 
fluctuations. This technique only measures the extent of 
clustered knots; stars or clusters in a diffuse component, 
even if originally born in clustered clusters, will evade de- 
tection. Therefore, we stress that any null detection does 
not indicate a lack of stellar populations at large galactic 
radii. 




Fig. 11. — The radii defined for creating the restricted three-point 
correlation maps. The center of the parent galaxy lies at the origin 
and the dot represents one outer disk knot. The origin al description 
of this method, along with this figure, are presented by Herbert-Fort 
|et al.|p009) . 

The restricted three-point correlation maps of our four 
low- inclination galaxies (i < 60°) are presented in Fig- 
ures [12] — [15] The top panels are constructed using all 
detections (the 'AH' sample) with -1.7 < U - V < 0.7 
and 19 < V < 27.5, while the lower panels result from 
splitting the samples into blue and red components on ei- 
ther side of U — V = —0.2 (middle and bottom panels, 
respectively, and comprise the 'Blue' and 'Red' samples). 
Black and gray show areas where signal is detected at the 
> 95% and > 99% significance level, respectively, as a 
function of both galactic radius R and intercluster radius 
Tout (Figure 11 ). The confidence l evels a re determined as 
described by |Herbert-Fort et al.| ( |2009[ ) using a Montc- 
Carlo approach. The dotted lines show the radial extent 
of H I for iV(H I) > 2 x 10 20 cm" 2 (left, at r out = kpc) 
and for N(R I) above the noise level of the integrated H I 
map (right, at r out = kpc). The dotted lines are slanted 
to distinguish the [R, r out ] regions that can be populated 
by sources within the H I disk. The H I data and analysis 
are described later in the paper. 



The interpretation of these Figures is somewhat unusual 
so we outline the salient features. Positive signal at any 
location is potentially a marker of clustering. However, 
positive signal at low r outl which is typically seen interior 
to i?25 is a sign of small-scale, tight clustering of stellar 
clusters. Vertical bands, which are seen at a variety of 
radii, suggest a set of clusters at a particular radius that 
may not be strongly clustered, such as those in a spiral 
arm or ring. Because of the nature of the fixed 
amount along the vertical axis represents a much larger 
angle at a smaller i?/i? 2 5 than at a larger one, but we 
expect clustering to depend on physical separation rather 
than on angular separation. 
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Fig. 12. — Restricted three-point correlation maps from sources 
in our (masked) final catalog of IC 4182. Black and gray show areas 
where signal is detected at the > 95% and > 99% significance level, 
respectively, as a fu ncti on of both galactic radius R and intercluster 
radius r ou t (Figure |ll|. The dotted lines show the radial extent of 



N(H I) above the noise level of the integrated H I map (right, at 
Tout = kpc) . The dotted lines distinguish the regions where knots 
could be co-located with the H I (left and right of the dotted lines, 
respectively). 
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Fig. 15.— Same as Figure \V2\ but for NGC 5474. 



Reaching conclusions regarding the radial extent of cor- 
related clusters from these maps is difficult because of the 
irregular nature of the signal and the statistical fact that 
given enough pixels some will randomly host 3<r variations. 
Rather than judge the significance of correlation peaks by 
eye, we calculate at each radius whether the occurrence of 
outliers is statistically beyond what is expected. For this 
calculation, we use probabilities calculated using binomial 
statistics with the appropriate probability corresponding 
to that of the unlikely event (for example, we calculate the 
likelihood of having as many as M or more pixels out of N 
pixels hosting 'true' values if the probability of obtaining 
a true value is 0.1). We rebin our data in the R direction 
so that each element now corresponds to a AR of the same 
size as the Ar ou t, so that each element is independent. 

We calculate the probability of finding as many 'signif- 
icant' detections in the correlation maps by chance as a 
function of galactocentric radius, over the range of r out 
plotted in the correlation maps (~ 5 kpc), and present 
the results in Figures 16 — 19. The black and gray lines 
correspond to the black and gray signal in the correlation 
maps. The dotted lines indicate the same H I extents as 
in the correlation maps. Black or gray signal below the 
dashed line at —2.568 indicates concentrations of correla- 
tion signal that are highly unlikely (> 3a) to be generated 
by chance (i.e. the correlation signal excess is statistically 
significant). 

Before discussing the results of our LBT correlation maps 
and the corresponding probability plots, wc present a sim- 
ilar correlation analysis of GALEX knots around these 
same galaxies and a cross-correlation analysis of our LBT 
knots and the underlying H I distribution. We then discuss 
the results of these three analyses together, on a galaxy-by- 
galaxy basis. We do this because comparison to indepen- 
dent data provides additional credence to our detections. 

4.2. Self- clustering of GALEX Knots 

We now apply our three-point correlation analysis to the 
distribution of UV-bright knots around our galaxies, using 
data from publicly available GALEX catalogs. If the outer 
disk regions that we detect in the LBT data are continu- 
ing to form clusters, then we might detect corresponding 
GALEX knot self-clustering (i.e. GALEX knot - GALEX 
knot clustering). Such a detection would not only con- 
firm what may be marginal detections in the LBT data, 
but also provide information on the timescales of spatially 
localized cluster formation and cluster dispersal. 

When possible, we combine the GALEX catalogs from 
shallow (t exp < 1000 seconds) and deep exposures to in- 
clude sources both near and far from the bright inner 
galactic disk (regions around the galaxies were masked 
in the catalog generation for the deep GALEX exposures 
and so the shallow exposures help fill in the source dis- 
tribution). We then cut the combined GALEX source 
lists t o match the 'blue' sample of Zaritsk y fc Christlein| 
( [2007] ) with FUV-NUV < 1 and NUV < 25, correspond- 
ing to clusters younger than ~ 360 Myr. This cut removes 
sources with detections in only one GALEX band and re- 
tains the most reliable knots. We reduce contamination 
further by only retaining knots that match sources in our 
LBT catalogs (to within 3", roughly half a GALEX resolu- 
tion element). The large majority, > 90%, of the GALEX 
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Fig. 16. — Probability, as a function of radius, tha t signal in 
the three-point correlation map of IC 4182 (Figure |12| is due to 
random excursions. The black and gray lines correspond to the black 
and gray signal in the correlation map. The dotted lines indicate 
the same H I extents as in the correlation map. Black or gray 
signal below the dashed li ne a t -2.568 indicates concentrations of 
correlation signal in Figure |l2"| that are very unlikely (> 3cr) to be 
generated by random noise. 
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Fig. 18.— Same as Figure 16 but for NGC 4736. 
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knots positionally match an LBT source. 

Before proceeding to the correlation analysis, we com- 
pare the number of GALEX knots around each galaxy 
with those remaining in the background-subtracted Hess 
diagrams (Figure 3 — 9) and present the results for two 
outer disk annuli in each low-inclination galaxy, between 
1.0 - 1.5R 25 and 1.5 - 2.0i? 25 , in Tabled Unfortunately, 
the inner annulus is problematic because the GALEX cata- 
logs near the disk are made from shallow exposures. In the 
1.5 — 2.0i?25 annulus, we find that GALEX knots make up 
between 10%-40% of LBT knots, or 27% on average. For 
a population of clusters that form continually, this ratio 
should simply reflect the timescales over which the clus- 
ters would be detected in GALEX versus LBT imaging, 
if cluster dissolution affects both samples similarly. Since 
dissolution is expe ct to proceed very rapidly (few m illion 
year timescale; see |Fall, Chandar, fc Whitmore|2005| ) both 
the GALEX and LBT samples, unlike Ha samples, are 
dominated by clusters that 'survive'. The average age of 
the GALEX population is ~ 180 Myr (assuming a uniform 
distribution of ages between — 360 Myr, where the upper 
bound is set by the UV color selection), while that of the 
LBT clusters is ~ 500 Myr (assuming a uniform distribu- 
tion between 0—1 Gyr, where the upper limit is again 
given by the color and magnitude selection). The ratio 
of these mean ages leads to an estimated fraction of 36%, 
sufficiently close to our detected average of 27% that we 
cannot rule out the hypothesis of continual cluster forma- 
tion. A lower fraction implies a higher cluster formation 
rate in the past, while a higher fraction implies a lower 
formation rate. The observed scatter between 10% - 40% 
suggests that factor of two variations in the cluster for- 
mation rate over these timescales are likely, but that the 
rates do not change by orders of magnitude over the pre- 
vious Gyr when integrated over timescales of ~ 200 Myr 
(roughly the resolution limit of our crude age estimates). 

We present the UV knot self-clustering maps from our 
final GALEX source catalogs in Figures [20] to [23j Unfor- 
tunately, the number of available GALEX sources is low; 
171 were used around IC4182, 179 around NGC 3351, 398 
around NGC 4736, and 174 around NGC 5474. Of these, 



NGC 4736 and NG C 5474 were included in the |Zaritsky 
& Christlein ( 2007 ) study. They found an excess of knots 
in both cases, although the probability of the excess be- 
ing random is not exceedingly small, 4% for NGC 4736, 
so the results based solely on the radial distribution of 
GALEX knots are marginal. Interestingly, the LBT re- 
sults for NGC 4736 are the strongest in our optical sample, 
and so show the value of going to deeper, redder samples 
in uncovering outer disk populations. 

4.3. Cross- correlation of LBT Knots and Neutral Gas 

Disks 

We now compare the distribution of our LBT sources 
to the neutral gas distribution. This cross-correlation pro- 
vides another way to constrain the extent of the stellar 
disks, particularly in relation to the known component 
that stretches to the farthest radii. Any association be- 
tween the H I structure and knots is also key for under- 
standing how star formation arises in this environment. 

Once again, because we perform an angular cross cor- 
relation analysis, we are restricted to our low-inclination 
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Fi g. 20. — Restricted three-point correlation map similar to Fig- 
ure [12] but here using 171 GALEX sources around IC 4182. The 
statistics are limited in comparison to the LBT self-clustering analy- 
sis by the relative shallower depth and poorer resolution of GALEX. 
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sample (IC 4182, NGC 3351, NGC 4736, and NGC 5474). 
A further limitation is that we do not have a suitable H I 
map of NGC 5474. The integrated H I (moment 0) maps 
for NGC 3351 and NGC 4 736 come from The H I Nearby 
Galaxy Survey (THINGS; [Walter et al.1[2008| and for IC 
4182 from the H I map of from the Wester bork H I Sur- 



vey of Sp iral and Irregular Galaxies (WHISP; Swaters et 



al. 2002). We estimate the extent of the gas in t he same 



manner as described by Herbert-Fort et al. (2009), by ex- 
amining histograms of the number of H I pixels per kpc 2 
above a certain N(H I) threshold, as a function of R in 
0.025i?25-wide elliptical annuli. 

The first threshold we consider is the noise level of the 
H I maps (see below) , to determine the maximum detected 
extent of the gas in each of the three galaxies. The 'noise 
level' of the final H I maps results from already having 
cleaned the maps, only keeping pixels where two or more 
adjacent channels show significant signal (2cr) in the inte- 
grated moment map (so the 'noise level' here is really 
more significant that the '2a' noise level in the raw data). 
The neutral gas around IC 4182 extends to <~ 2.25i?25 
around IC 4182, to - 2.4i? 25 around NGC 3351, and to 
- I.8-R25 around NGC 4736. Figures [24] — [26] show the 
restricted three-point cross-correlation between the LBT 
knots and the H I pixels lying above the A(H I) threshold 
(listed in the legend at the bottom right of each panel). 
The top panels show results using the 'All' sample while 
the lower panels give the results from using the 'Blue' and 
'Red' samples (middle and bottom panel, respectively). 
Black and gray show areas where signal is detected at 
> 95% and > 99% significance, respectively. The dotted 
lines bracket the radial extent of the H I pixels used. The 
Figures show that LBT knots correlate with H I features 
ou t to, and beyond, the obs erved edges of the gas disks. 
As Herbert-Fort et al. ( 2009[ ) discussed, the vertical bands 
of signal highlight likely spiral arm structures. 

To determine what H I density might be the best tracer 
of knot formation, we focus on our 'Blue' LBT knot sam- 
ple and calculate correlations relative to the H I distri- 
bution as defined by different column density thresholds. 



Figures 27 — 29 contain results specifically for AT(H I) > 



1.0 x lO^cm "middle), and A(H I) > 2.0 x 10 2 
(bottom). The highest threshold, N(H I) > 2.0 x 10 20 cm- 
corresponds to where one finds damped Lyman-a (DLA) 
absorption, and so, distinguishes region s that contain pre- 



dominantly neutral or ionized gas (see Wolfe et al.| 2005 
and references therein). The A(H I) = 2 x 10^ u cm"' 



1). The N(H I) 
threshold highlights the edge of the dominant reservoir 
of neutral gas. The N(R I) = 1.0 x 10 20 cm" 2 threshold 
was determined by eye to be the lowest density value at 
which the contours closely trace the distribution of the 
'Blue' knots. The radial extent of the H I decreases as the 
column density threshold increases, though as the den- 
sity threshold is increased we are more likely selecting H I 
structures that could host cluster formation — unless the 
gas has been consumed to make the clusters or dispersed 
after cluster formation. Broadly, we find correlations be- 
tween the knots and the gas at all gas densities and at all 
radii. We discuss individual cases and reach conclusions 
next. 



5. RESULTS FOR INDIVIDUAL GALAXIES 
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Fig. 24. — Restricted three-point cross-correlation maps of LBT- 
detected knots around IC 4182 and H I pixels with 7V(H I) above 
the noise level of the integrated H I map (iV(H I) > 4.0 X 10 16 cm -2 
here). The dotted lines bracket the radial extent of the H I pixels 
used; because no H I pixels exist beyond the dotted line on the right, 
we will not see signal at low r ou t beyond this furthest H I radius 
(this explains the diagonal nature of signal at the largest R). Only 
knots between 1.0 — 3.O-R25 were used. 
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Fig. 25. — Same as Figure [24| but for NGC 3351. Only sources 
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A few general comments before we discuss the galaxies 
individually. In all four LBT correlation maps (Figures 12 



15) there is significant signal at very large radii, in 



three cases extending beyond the H I (we do not have H I 
data for the fourth galaxy, NGC 5474). Though we are 
somewhat skeptical of correlation signal at such large radii, 
where the background becomes relatively more important, 
these detections may point to stars formed beyond the cur- 
rent edge of the observable H I disk. We examined the 
galaxy images and, with one exception that we discuss be- 
low, could not visually identify groupings of LBT knots 
at radii beyond the H I extent. However, the knots can 
be very difficult to find when they are far from the disk, 
not strongly clustered, and faint, which is of course why 
we devised the restricted 3-pt correlation analysis. When 
we compare correlation maps made using the masked and 
unmasked source distributions (masked to limit the cata- 
log problems near bright stars, as mentioned in §2), much 
of the signal beyond the H I changes significantly in char- 
acter, suggesting it might be artificial. We hesitate to in- 
crease the masking beyond what we already chose because 
of the limitation in the number of sources near the mosaic 
edges. Nevertheless, although some of the most distant 
signal may be artificial, we will demonstrate that some of 
the features are real and that there are clusters out to at 
least ~ 3i?25 and beyond the edge of the H I disk. We now 
discuss the results for each galaxy. 

5.0.1. IC 4182 

We find the bulk of the signal contained within ~ R25 in 
the LBT knot correlation map, although there is a great 
deal of structure in the Blue knot correlation map at larger 
radii, and some conc ent ration of signal toward low r out in 
the Red one (Figure 12). Working our way out in radius, 
the first tantalizing detection in the 2-D plot is in the Red 
knots at about ~ 1.5i? 2 5- Howeve r th is turns out to not 
be statistically significant (Figure 4.2) and is also likely 



to be the remnant of structure caused by a bright star 
at that radius. However, strong signal exists near 1.5i?25 
in the All and Blue panels of the LBT knot -HI cross- 
correlation plot (Figure 24 ) . The lack of a self-clustering 
signal among the Blue sample is not counter-evidence, just 
evidence for a lack of strong self-clustering among such 
clusters. Interestingly, the Blue LBT knot -HI cross 
correlation at 1.5i?25 is most closely related to the lower- 
density gas, N(H I) < 20.0 dex (Figure 27), suggesting 
that these knots have either moved away from the regions 
of highest gas density or affected their nearby environs. 

Further out in R, the Blue panel of the LBT knot prob- 
ability plot (Figure 16) shows more noisy, scattered signal 
beyond the edge of the H I than do the All or Red panels. 
This 'noisy' Blue signal in the LBT correlation map be- 
tween 2 — 3i?25 is statistically significant in the probability 
plot and coincides with bands of Blue signal in the LBT 
knot -HI cross-correlation plots, suggesting that Blue 
clusters a few kpc beyond the edge of the gas are neverthe- 
less strongly associated with the outer disk gas structure. 
This result suggests that the knots were either born from 
gas that was once part of the disk, but no longer exists, or 
that they have drifted somewhat in radius from their birth- 
sites. A drift velocity of just 6 km s _1 could transport a 
cluster a few kpc in 500 Myr. The GALEX correlation plot 



(Figure 20) does not contain signal at these radii, which 
suggests that these clusters are somewhat older than a few 
hundred Myr. This scenario is consistent with the infer- 
ences discussed above that the clusters are not as strongly 
correlated with the highest density gas and that they may 
have drifted in radius. 

Finally, all of the panels in the LBT knot correlation 
probability plot (Figure 16) show significant excess signal 
at ~ 3.3i?25, far outside the edge of the gas disk (the gas 
extends to 2.2i? 2 s)- This feature is especially interesting 
because it corresponds to strong signal at the sa me r adius 
in the GALEX knot self-correlation map (Figure 20). Ex- 



amining the GALEX image, we visually identify a clump 
of sources at that radius. If these distant sources are in 
fact outer-disk clusters, rather than a grouping of UV- 
bright background objects or stellar clusters in a satellite 
galaxy, they are very interesting objects for further study. 
Ha imaging, and subsequent spectroscopy would be par- 
ticularly valuable. 

5.0.2. NGC 3351 

The bulk of the signal in the LBT knot corr elat ion map 
and corresponding probability plot (Figures 13 and 17) 
is once again contained within ~ R25, demonstrating that 
the optical radius does demarcate a real qualitative change 
in the nature of star formation. Again moving out in ra- 
dius, the most significant signal comes in at around 2.4i? 2 5 
in the All and Red panels, and to a lesser degree in the 
Blue panel. These clusters are also seen in the GALEX 
correlations (Figure 21). Evidence for correlation between 
this population and the H I is weaker, although there is 
some signal near 2.2_R 2 5 i n the Blue panel and 2.5i?25 in 



the Red panel (Figure 25 ) . This signal occurs right at the 
edge of the H I disk. 

In addition to the clusters seen at the periphery of the 
H I, there is also signal in the LBT knot -HI cross- 
correlation map form R25 out to I.6JZ25. From the Blue 
LBT knot -HI cross-correlation map for different iV(H I) 
thresholds (Figure 28), we conclude that the knots trace 
the denser gas out to ~ I.6A25 (and trace spiral structure 
out to 1.2i?25)- As in IC 4182, the bands of correlations 
between the Blue knots and H I come at semi-regular ra- 
dial spacings, suggestive of multi-armed spirals with small 
pitch angles or rings. 

5.0.3. NGC 4736 

As we have discussed before, NGC 4736 is by far the 
richest of our galaxies in outer disk structure and may 
therefore not be representative. Nevertheless, we see many 
of the same general patterns. The LBT knot correlation 
map and corresponding probability plot for NGC 4736 
(Figures 14 and 18) show signal extending to at near the 
limit of the HI disk (~ I.6.R25). This galaxy has an ob- 
vio us overdensity of kno ts in the bright 'ring' surrounding 
it. |Trujillo et al. (2009) found that this ring is actually a 



complex structure of wound spiral arms (perhaps caused 
by secular processes in the inner disk), which perhaps are 
a much stronger version of the low pitch angle arms we 
discussed in the context of IC 4182 and NGC 3351. The 
knots toward the inner edge of the outer disk are associ- 
ated with higher density H I, as is also the case for NGC 
3351. We are somewhat skeptical of excess signal beyond 
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the H I extent of this galaxy (signal at R ~ 2i?2s) be- 
cause it is coincident with artificial signal that faded when 
we applied additional bright star masking, though there is 
some signal in the H I correlation maps at the same ra- 
dius. The lack of correlation in the GALEX data is not 
a concern because the signal in the LBT correlation maps 
comes primarily from the Red sample. 

5.0.4. NGC 5474 

NGC 5474 is our lowest noise LBT knot correlation map 
and contains significant signal from R25 to ~ IAR25 (Fig- 
ures 15 and 19). Unfortunately, we do not have the same 
quality H I data for NGC 5474 as for the other three galax- 
ies. However, the GALEX knot correlation map contains 
some signal near 1.5i?25, though it is very weak (Figure 
23 ) . There are two reasons why the signal might be weak. 
First, it is possible that the GALEX knots do not cluster 
very tightly (or do not cluster differently than the back- 
ground). Second, it may be that the relative small number 
of GALEX knots limits the degree to which this restricted 
3-pt correlation can result in a statistically meaningful de- 
tection. From looking at the GALEX images and source 
distribution we conclude that this feature is real, but this 
weak signal (as well as the rather tame GALEX signal for 
NGC 4736) demonstrates how our optical work comple- 
ments the GALEX studies. 

There are a few marginally significant detections at in- 
termediate radii in the correlation probability plot (Figure 
19). A potentially interesting detection is that in the All 
panel at a radius slightly below 2i?25- This radius matches 
our est imate for the outer extent of H I based on Figure 



2 from|Rownd et al. (TL994|). While the other three galax- 
shc 



ies in our sample all show clear signs for cluster formation 
near the periphery of their H I distribution, we only have 
this questionable detection for NGC 5474. 

In the outer extremity of the correlation maps, R ~ 3.5 
to 4 R25, we again find strong signal in the three knot 
populations. There is no corresponding signature in the 
GALEX maps (although we saw above that the GALEX 
detections can be weak) and, unfortunately, we do not 
have H I maps of the required quality to help confirm 
these features. The feature near 4i?25 is likely residual 
signal from a few bright stars at that radius. 

6. SUMMARY AND CONCLUSIONS 

We summarize our results as follows: 

• All disk galaxies have a cluster population beyond 
i?25- The six for which wc present background subtracted 
Hess diagrams show significant populations out to at least 
1.5i?2 5- The outer disk we stud ied similarly before (NGC 
3184; |Herbert-Fort et al. 2009) also contains this popula- 
tion. Many of our galaxies show similar populations ex- 
tending to 2i?25 and occasionally beyond. We attribute the 
larger fraction of galaxies with detec ted outer disk popu- 



lations relative to previous studies (Thilker et al. 



2007 
d age 



Zaritsky & Christlein 20071 to the superior mass an 
sensitivity of the LBT data. 

• Using the distribution of sources in the Hess diagrams, 
we infer that the typical detected cluster has a mass of 
~ 1O 3 M , is predominantly < 1 Gyr old, and as a pop- 
ulation have an average formation rate of at least ~ one 
cluster every 2.5 Myr. The corresponding rate of stellar 



mass being formed in these clusters is ~ 0.004M Q pc -2 
Gyr -1 (assuming R25 = 5 kpc and 10 3 M Q clusters) for 
the area of the disk between R25 and 1.5i?25- These es- 
timates are rough, and may be systematically biased due 
to uncertainties in the modeling, stochastic effects in low 
mass cluster, and selection. The principal aim of this exer- 
cise was to demonstrate the plausibility of associating the 
detections with stellar clusters. Comparing the numbers 
of sources identified to the corresponding numbers found 
using GALEX, which is sensitive only to younger clusters, 
we conclude that the formation rate of clusters at these 
radii is constant, to within a factor of roughly 2, over the 
last Gyr. 

• To further quantify the distribution of outer disk clus- 
ters, we construct restricted three-point correlation maps 
in our four low-inclination galaxies. We confirm many of 
these detections using comparisons of GALEX detected 
populations and correlations with H I. Again, we detect 
signal in all four, but our detections come in three vari- 
eties. First, we generally find a population of clusters that 
extends modestly beyond the optical radii (to between 1.3 
and 1.5 -fi^s)- Second, we find a population of clusters 
near the edge of the H I distribution. Lastly, in all but 
NGC 3351, we find detections of clusters well beyond the 
H I edge. These last are the most difficult to confirm inde- 
pendently (they could either be an unfortunate clustering 
of background sources or they could belong to satellite 
galaxies) . 

• From the cross-correlation signal between our Blue 
LBT-detected knots and the H I distribution we find two 
types of behavior. First, the knots near the optical edge of 
the disk are best traced by the higher density neutral gas, 
as one might expect in a steady state situation where den- 
sity waves continually lead to the generation of new clus- 
ters. The pattern of the correlations, semi-regular bands 
in Figures [24] to [26] also suggest spiral arms, in this case 
with small pitch angles. Second, the knots farther out in 
the disk are most strongly correlated to the low density H I 
gas, suggesting that while some connection between star 
formation and fuel exists, the process is sufficiently tran- 
sient and/or disturbing that correlations with high density 
gas do not persist. 

From these results, as well as those presented in previ- 
ous studies of outer disks, we suggest that outer disk clus- 
ter formation occurs in three modes. First, spiral waves 
from the inner disk continue beyond the optical radius and 
trigger continual, but low level, cluster formation out to at 
most ~ 1.5i?25- We have presented some evidence for this 
mode of cluster formation, but it can also be clearly seen in 
the Ha images of NGC 628 (Ferguson et al. 1998) where 
the arms can be visually traced beyond R2$- Second, a 
global mode, where clusters are formed throughout the 
disk, is triggered by interactions. This mode is responsi- 
ble for the rare and most dramatic example s of outer disk 
star f ormation, such as that seen in M83 (Thilker et al. 
2005). Lastly, and most speculative, is a mode that cre- 



ates clusters at the periphery of the H I disk. We suggest 
that this is where gas that is accreted joins the already 
existing gas and that this process leads to low level clus- 
ter formation. We consider these radii to represent the 
outer banks of galaxies, where the waves of incoming gas 
break upon the shores of the existing disk and leave telltale 



'foam' upon their arrival. 

This work relied on the invaluable idlutils software li- 
brary developed by M. R. Blanton, S. Buries, D. P. Finkbeiner, 
D. W. Hogg, and D. J. Schlegel, and on the Goddard 
IDL library maintained by W. Landsman. We thank the 
THINGS and WHISP teams that produced the H I maps 
used here. We thank the referee for suggestions that clari- 
fied the text. DZ and SHF were partially supported under 
NASA LTSA NNG05GE82G and NSF AST-0307482. DZ 
thanks Cambridge University and New York University for 
their hospitality during the final stages of this work. 



Table 1 



Sample Galaxies 



Name 




PA 


D 


R25 


Udate 


Uexp 


Vdate 




Vlim 


Vacor 


a v acoT 


Ref. 




[deg] 


[deg] 


[Mpc] 


[kpc] 




[sec] 




[sec] 








IC 4182 


23 


90 


4.7 


4.1 


5/11/07 


1640 


5/11/07 


1640 


27.0 


-0.11 


0.05 


1 


NGC 3351 


10 


13 


10.1 


10.8 


3/21/07 


1640 


3/21/07 


1640 


27.0 


-0.15 


0.05 


4,5,6,7 


NGC 4736 


8 


122 


4.7 


7.6 


2/21/07 


1640 


2/21/07 


1476 


27.5 


-0.13 


0.03 


2,3 


NGC 4826 


61 


115 


7.5 


10.2 


2/10/08 


1476 


4/24/07 


1476 


26.0 


-0.69 


0.04 


1,8,9 


NGC 5474 


26 


132 


7.2 


5.0 


5/10/07 


328 


5/10/07 


1804 


26.0 


-0.26 


0.01 


1,10 


NGC 6503 


74 


121 


5.3 


4.5 


4/23/07 


1312 


4/11,23/07 


3280 


26.0 


-0.46 


0.01 


1,11,12 



References. — 1: Karachcntscv ct al. (2004); 2: Karachcntscv (2005); 3: Tmjillo ct al. (2009); 4: Graham ct al. (1997); 5: Rubin 



et al.| 11975 


i; 6: |ButaT 


1988||; 7: |Kwartz et al.| 


Makarova (1999); 12: E 


Segeman et al.|(|1987p 



(2006); 8: Tonry ct al. (2001); 9: [ 



Nilson ct al. 



10: Rownd ct al. (1994) 



ubin| 



Table 2 

Sources in Color-Magnitude and Background-Subtracted Hess Diagrams 



Name 


Annulus 


Area 


NcMD 




Fraction 


Me// 






IBM] 


[arcmin - 2 ] 






[%] 


[mag arcsec ] 


IC 4182 


1.0 


- 1.5 


30.95 


1950 


248 


13 


29.62 




1.5 


- 2.0 


41.66 


2244 


-48 




30.44 




2.0 


- 2.5 


43.91 


2408 


-7 




29.98 


NGC 3351 


1.0 


- 1.5 


40.22 


2037 


448 


22 


> 32.44 




1.5 


- 2.0 


38.22 


1593 


8.3 


5 


> 34.91 




2.0 


- 2.5 


37.52 


1554 


72 


5 


> 31.12 


NGC 4736 


1.0 


- 1.5 


85.38 


8517 


2774 


33 


29.46 




1.5 


- 2.0 


104.88 


7570 


516 


7 


30.43 




2.0 


- 2.4 


83.93 


6135 


489 


8 


> 30.99 


NGC 4826 


1.0 


- 1.5 


40.42 


1445 


-33 




30.00 




1.5 


- 2.0 


52.90 


1964 


30 


2 


33.11 




2.0 


- 2.5 


72.15 


2670 


33 


1 


31.34 


NGC 5474 


1.0 


- 1.5 


19.96 


699 


427 


61 


28.85 




1.5 


- 2.0 


27.95 


433 


52 


12 


> 33.57 




2.0 


- 2.5 


35.93 


465 


-25 




32.14 


NGC 6503 


1.0 


- 1.5 


8.69 


622 


369 


59 


28.36 




1.5 


- 2.0 


9.94 


554 


265 


48 


28.94 




2.0 


- 2.5 


13.63 


511 


114 


22 


29.75 



Fraction is undefined when N[j eaa is negative (there is no excess of disk sources over 
the background). £t e // is given as a lower limit when the total flux in the Hess diagram is 
negative. Note that there exist cases where Nj{ esa is negative but the total flux is not. This 
depends on the distribution of signal across the entire Hess diagram. 



Table 3 

Sources in LBT-matched GALEX catalogs 



Name 


Annulus 


Ngalex 




Fraction 












[%] 


IC 4182 


1.0 


- 1.5 


5 


348 


1 


IC 4182 


1.5 


- 2.0 


28 


69 


41 


NGC 3351 


1.0 


- 1.5 


12 


548 


2 


NGC 3351 


1.5 


- 2.0 


27 


266 


10 


NGC 4736 


1.0 


- 1.5 


56 


2738 


2 


NGC 4736 


1.5 


- 2.0 


172 


704 


24 


NGC 5474 


1.0 


- 1.5 


8 


427 


2 


NGC 5474 


1.5 


- 2.0 


16 


52 


.31 



Because the GALEX catalogs already had so few usable 
sources we did not mask the regions we did in our LBT mosaics 
around bright stars; therefore, so that the fields are of similar 
area, the Nn es3 listed here are the values from our LBT cata- 
logs before the bright star masks were applied. 
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